Tunable "Doniach Phase Diagram" for strongly-correlated nanoclusters 
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Exact diagonalization calculations reveal that the energy spacing A in the conduction band tunes 
the interplay between the local Kondo and non local RKKY interactions, giving rise to a "Doniach 
phase diagram" for a nanocluster with regions of prevailing Kondo or RKKY correlations. The parity 
of the total number of electrons alters the competition between the Kondo and RKKY correlations. 
This interplay may be relevant to experimental realizations of small rings or quantum dots with 
tunable magnetic properties. Below a critical value V c of the hybridization the susceptibility exhibits 
a low-T exponential activation behavior determined by the interplay of the spin gap and A. 



The interplay between the local Kondo interactions and 
the non local Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interactions is a central unresolved problem in the 
physics of dense valence fluctuation and heavy fermion 
compoundspi*Sii The first interaction, responsible for the 
quenching of the local /-moment via the screening of the 
conduction electrons, has a characteristic energy scale 
given by the single-impurity bulk Kondo temperature, 
T K oc exp{—l/p{Ep)J)A Here, p(E F ) is the density of 
states of the conduction electrons at the Fermi energy 
and J is the antiferromagnetic interaction between the 
impurity spin and the conduction electrons. The latter 
interaction is an indirect magnetic interaction between 
localized moments mediated by the conduction electrons, 
with an energy scale of order Tukky oc J 2 p(Ep) which 
promotes long- or short-range magnetic orderingi£i4 In 
the high-temperature regime the localized moments and 
the conduction electrons retain their identities and inter- 
act weakly with each other. At low-temperatures, the 
moments order magnetically if the RKKY interaction is 
much larger than the Kondo energy, while in the reverse 
case the system forms a heavy Fermi liquid of quasipar- 
ticles which are composites of local moment spins bound 
to the conduction electrons^ The overall physics can be 
described by the well-known "Doniach phase diagram" i 
The description of the low-temperature state, when both 
the RKKY and the Kondo interactions are of comparable 
magnitude, is an intriguing question that remains poorly 
understood and is the subject of current active research. 3 

There has been a revival of the Kondo physics in 
recent years due to an exciting series of experiments 
on nanoscale systems, which has enabled to probe 
the local interactions in a well controlled way at the 
nanoscalei£i£iL£ For example, STM experiments^ stud- 
ied the interaction of magnetic impurities with the elec- 
trons of a single-wall nanotube confined in one dimen- 
sion. Interestingly, in addition to the bulk Kondo res- 
onance new sub peaks were found in shortened carbon 
nanotubes, separated by about the average energy spac- 
ing, A, in the nanotube. These experiments invite an 
interesting and important question: How is the inter- 
play between the Kondo effect and RKKY interactions 
becomes modified in a nanocluster containing magnetic 
impurities, where the conduction electron spectrum be- 
comes discrete with a mean energy level spacing A. 



For such small systems, controlling Tk upon varying 
A is acquiring increasing importance since it allows to 
tune the cluster magnetic behavior and to encode quan- 
tum information. While the effect of A on the single- 
impurity Anderson or Kondo model has received consid- 
erable theoretical2ii iii*i2ii 3 rii4attention recently, its role 
on dense impurity clusters remains an unexplored area 
thus far. The low-temperature behavior of a nanosizcd 
heavy-electron system was recently studied within the 
mean-field approximation^. 

In this work we present exact diagonalization calcula- 
tions for d- or /-electron nanoclusters to study the ef- 
fect of energy spacing, parity of number of electrons, and 
hybridization on the interplay between the Kondo and 
RKKY interactions in dense strongly correlated electron 
nanoclusters. While the properties of the system depend 
on their geometry and sizoi&, the present calculations 
treat exactly the Kondo and RKKY interactions. Our 
results show that: i) tuning A and the parity of the total 
number of electrons can drive the nanocluster from the 
Kondo to the RKKY regime, i.e. a zero- temperature en- 
ergy spacing versus hybridization phase diagram; ii) the 
temperature versus hybridization "Doniach" phase dia- 
gram for the nanocluster depends on the energy spacing 
; and iii) Below a critical value V c of the hybridization 
the susceptibility exhibits a low-T exponential activation 
behavior determined by the interplay of the spin gap and 
A. In contrast, above V c there is no exponential behav- 
ior. 

The one dimensional Anderson lattice Hamiltonian is 
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where, t is the nearest-neighbor hopping matrix element 
for the conduction electrons (e/. = —2tcosk), (cj )CT ) 
annihilates a localized / (conduction) electron on site i 
with spin a, Ef is the energy level of the localized orbital, 
V is the on-site / — c hybridization matrix element and 
U is the on-site Coulomb repulsion of the f electrons. 
We consider the half-filled (N e i = 2N = 12) symmetric 
(Ef = — y = —2.5) case. As expected, the cluster has a 
singlet ground state (S g = where S g is the ground-state 
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spin) at half filling. The energy spacing in the conduction 
band is A = At/ (N — 1 ) = 4p The exact diagonalization 
calculations employ periodic boundary conditions. We 
compare the onsite local Kondo spin correlation function 
(SCF) < S{ z Sf z > and the nonlocal (nearest-neighbor) 
RKKY SCF < Sl z S f (i+1) z > to assign a state to the 
Kondo or RKKY regimes, in analogy with mean field 
treatments H 

In Fig. 1 we present the variation of the local Kondo 
SCF < S{ Z S? Z > (squares) and the nearest-neighbor 
RKKY SCF < S[ z S f (i+1) z > (circles) as a function of 
hybridization for two values of the hopping matrix ele- 
ment t = 0.2 (closed symbols) and t = 1.2 (open sym- 
bols), respectively. For weak hybridization the nonlo- 
cal (local) nearest-neighbor RKKY (Kondo) SCF is large 
(small), indicating strong short-range antiferromagnetic 
coupling between the f — f local moments, which leads 
to long range magnetic ordering for extended systems. 
As V increases, < S{ tZ St +1 ^ z > decreases whereas the 

< S( z Sf z > increases (in absolute value) saturating at 
large values of V. This gives rise to the condensation of in- 
dependent local Kondo singlets at low temperatures, i.e., 
a disordered spin liquid phase. For large V the physics 
are local. Interestingly, as t or A decreases the f-c SCF is 
dramatically enhanced while the /-/SCF becomes weaker, 
indicating a transition from the RKKY to the Kondo 
regime. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 



FIG. 1: Nearest-neighbor f-f SCF (circles) and on-site f-c SCF 
(squares) as a function of V for t = 0.2 (closed symbols) and 
t — 1.2 (open symbols), respectively. 

In Fig. 2 we present the energy spacing versus V zero- 
temperature phase diagram of the nanocluster, which il- 
lustrates the interplay between Kondo and RKKY in- 
teractions. In the RKKY region < S{ z S^ i+1 ^ z > is 

larger than < S{ z Sf z > and the total local moment 
squared [i 2 = ((^/ + f-c) 2 ) 0; m the Kondo regime 



< S , { z S , ( / i+1) z > is smaller than the < Sl z Sf z >, fi 2 = 0, 
and the ground state is composed of independent local 
singlets. The solid curve indicates the critical value of 
the hybridization V — V c or the critical value of energy 
spacing A = A c , where the local and non local SCF's are 
equal, i.e., < S f lz S f {i+1)z >=< S{ z Sf z >. The dashed 
curve denotes the points where the total local moment 
square vanishes. Thus, in the intermediate regime, the 
so-called free spins regime^, < S{ z S^ i+l ^ z > is smaller 

than the < S( z Sf z > , the / moment is partially quenched 
and [x 2 ^ 0. Interestingly, we find that the free spins 
regime becomes narrower as the average level spacing A 
is reduced. This result may be interpreted as a quan- 
tum critical regime for the nanocluster due to the finite 
energy spacing, which eventually reduces to a quantum 
critical point when A — > 0. 
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FIG. 2: Energy spacing A versus hybridization zero- 
temperature phase diagram. The solid curve denotes the 
crossover point of the spin-spin correlation function in Fig.l; 
the dashed curve denotes the set of points where the on-site 
total moment square (([if + u c ) 2 ) = 0.0 ± 0.05. 

The T=0 exact diagonalization results on small clus- 
ters are generally plagued by strong finite size effect s^SiiS. 
Performing calculations at finite T gives not only the 
thermodynamic properties of the system, but most im- 
portantly diminishes finite-size effects for ksT ^> A. In 
Fig. 3 we show the nearest-neighbor f-f SCF and on- 
site Kondo f-c SCF as a function of temperature for 
t = 0.2 and V = 0.2 < V c and V = 0.4 > V c , where 
V c = 0.25. At high temperatures, the free moments of 
the f and conduction electrons are essentially decoupled. 
The nearest-neighbor non local SCF falls more rapidly 
with T than the on-site local f — c SCF, indicating that 
the non local spin correlations are destroyed easier by 
thermal fluctuations. For V < V c , the nanocluster is 
dominated by RKKY interactions at temperatures lower 
than the crossover temperature, T^ KKY , which denotes 
the temperature where the non local and local interaction 
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become equal in the nanocluster. In the infinite system 
this temperature would denote the ordering Neel temper- 
ature. On the other hand, for V > V c the RKKY and 
Kondo spin correlation functions do not intersect at any 
T, and the physics are dominated by the local Kondo 
interactions. 




100 



T(K) 



FIG. 3: Nearest-neighbor f-f and on-site f-c SCF versus tem- 
perature for t = 0.2 and V = 0.2 < V c and V = 0.4 > V c , 
V c = 0.25. 
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FIG. 4: Effect of energy spacing, A = t^t on the exact 
"Doniach phase diagram" for a strongly correlated electron 
nanochain. The curve represents the crossover temperature 
Trkky i where the non local short range AF spin correlations 
become equal to the local on-site Kondo spin correlations. 
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In Fig. 4 we present the crossover temperature T^ KKY 
for the cluster as a function of hybridization for different 
values of t. This represents the phase diagram of the 
strongly correlated nanocluster, which is similar to the 
"Doniach phase diagram" for the infinite Kondo necklace 
modeli The phase for T < T^ KKY denotes the regime 
where the non local short-range magnetic correlations are 
dominant. For V < V c and T » T^ KKY one enters 
into the disordered "free" local moment regime. On the 
other hand, for V > V c and at low T, the nanocluster 
can be viewed as a condensate of singlets, typical of the 
Kondo spin-liquid regime. Interestingly, the T^ KKY can 
be tuned by the energy spacing A or the size of the clus- 
ter. Thus, increasing A or decreasing the size of the nan- 
ocluster results to enhancement of the non local nearest- 
neighbor RKKY correlations and hence Tfj KKY . This 
result is the first exact "Doniach phase diagram" for a 
nanocluster. 

In bulk Kondo insulators and heavy- fermion systems, 
the low-T susceptibility and specific heat behavior is de- 
termined by the spin gap, which for the half-filled An- 
derson lattice model, is determined by the ratio of V to 
U. On the other hand, strongly correlated nanoclusters 
are inherently associated with a new low-energy cutoff, 
namely the energy spacing A of the conduction electrons. 
Thus, a key question is how can the low-temperature 
physics be tuned by the interplay of the spin gap and the 
energy spacing. The temperature-dependent local f-spin 



where Q is the partition function. Here, S is the 
z-projection of the total spin (both the /- and o 
contributions), and \a > and E a are the exact many- 
body eigenstates and eigenvalues, respectively. When 
V = 0, the localized spins and conduction electrons are 
decoupled and X ^(T) is simply the sum of the Curie 
term due to the free f spins and the Pauli term of the 
free conduction electrons. For finite V, X CO decreases 
with temperature at low-temperatures. In Fig. 5 we 
present x CO as a function of temperature for t = 0.2 
and V = 0.2 < V c , V = V c = 0.25, and V = 0.4 > V c . 
For small V the spin gap which is smaller than A controls 
the exponential activation behavior of \ &t l° w ^ ( see 
inset in Fig. 5). In contrast, in the large V limit, the 
spin gap becomes larger than A and the low-T behav- 
ior of the susceptibility shows no exponential activation. 
At high T we can see an asymptotic Curie- Weiss regime, 
typical of localized decoupled moments. The low-T be- 
havior of the susceptibility is associated with the lowest 
energy scale determined by the lowest value between the 
spin gap and the energy spacing A. For large values of 
t (or A) the spin gap is reduced and the spin gap is the 
lowest energy scale. Consequently, the low-T behavior 
exhibits exponential activation associated with the spin 
gap. Similar low-T behavior is also found for the specific 
heatiSi 

We have also examined the role of even versus odd 
number of electrons on the magnetic behavior. For t = 1, 
changing the number of electrons from N e ; = 12 to N e ; 
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FIG. 5: Local f magnetic susceptibility as a function of tem- 
perature for t = 0.2 and V = 0.2 < V c , V = V c = 0.25 and 
V — 0.4 > V c . The inset shows the linear fit for ln(x ) versus 
i at the low temperature. 

= 11 results in: (a) an enhancement of the local Kondo 
f-c spin correlation function; and (b) a suppression of the 
nearest-neighbor f-f spin correlation function. This in- 
teresting novel tuning of the magnetic behavior can be 
understood as follows: For an odd-electron cluster, the 
topmost occupied single particle energy level is singly 
occupied. On the other hand, for an even-electron clus- 
ter, the topmost occupied single-particle energy level is 
doubly occupied, thus blocking energy-lowering spin-flip 



transitions. This energy penalty intrinsically weakens the 
Kondo correlations^. As expected, changing the number 
of electrons from even to odd changes S g = to S g = h. 
For t= 0.05 (Kondo regime), the on site f-c correlation 
function does not depend as strongly on the parity in the 
number of electrons because the sites are locked into sin- 
glets. On the other hand, the nonlocal RKKY becomes 
suppressed as in the case of large energy spacing. 

Exact diagonalization calculations reveal that the en- 
ergy spacing and the parity of the number of electrons, 
give rise to a novel tuning of the magnetic behavior of 
a dense strongly correlated nanocluster. This interesting 
and important tuning can drive the nanocluster from the 
Kondo to the RKKY regime, i.e. a tunable "Doniach" 
phase diagram in nanoclusters. For weak hybridization, 
where the spin gap is smaller than A, both the low- 
temperature local / susceptibility and specific heat ex- 
hibit an exponential activation behavior associated with 
the spin gap. In contrast in the large hybridization limit, 
A is smaller than the spin gap, the physics become local 
and the exponential activation behavior disappears. We 
believe that the conclusions of our calculations should be 
relevant to experimental realizations^ of small clusters 
and quantum dots. 
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